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HIGH ENERGY NEUTRINOS AND COSMIC-RAYS 
FROM LOW-LUMINOSITY GAMMA-RAY BURSTS? 

KOHTA MURASE 1 , KUNIHITO IOKA 2 , SHIGEHIRO NAGATAKi'- 3 , AND TAKASHI NAKAMURA 2 
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ABSTRACT 

The recently discovered gamma-ray burst (GRB) GRB 060218/SN 2006aj is classified as an X-ray Flash with 
very long duration driven possibly by a neutron star. Since GRB 060218 is very near ~ 140 Mpc and very dim, 
one-year observation by Swift suggests that the rate of GRB 060218-like events might be very high so that such 
low luminosity GRBs (LL-GRBs) might form a different population from the cosmological high luminosity 
GRBs (HL-GRBs). We found that the high energy neutrino background from LL-GRBs could be compara- 
ble with that from HL-GRBs. If each neutrino event is detected by IceCube, later optical-infrared follow-up 
observations such as by Subaru and HST have possibilities to identify a Type Ibc supernova associated with 
LL-GRBs, even if gamma- and X-rays are not observed by Swift. This is in a sense a new window from neutrino 
astronomy, which might enable us to confirm the existence of LL-GRBs and to obtain information about their 
rate and origin. We also argue LL-GRBs as high energy gamma-ray and cosmic-ray sources. 

Subject headings: gamma rays: bursts — acceleration of particles — elementary particles 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) and supernovae (SNe) are most 
powerful phenomena in the universe. Theorists predicted that 
the former would result from the death of massive stars, and 
the association of long-duration GRBs with core-collapse su- 
pernovae (SNe of Type Ibc to be more specific) has been ob- 
served over the last decade. The first hint for such a connec- 
tion came with the discovery of a near by SN 1998bw (SN 
Ic) in the error circl e of GRB 980425 dGalama et al. 19981 
Iwam otoetal. 199 81 The first spectroscopic identification of 
a SN Ic superposed on a GRB after glow compon ent was done 
in GRB 030329/SN 2003dh (Hiort h et al 20031) . 

Recently Swift discovered GRB 060218, which is the sec- 
ond nearest GRB identified to-date (Campana et al. 2006; 
Cusumano et al. 2006; Mirabal & Halpern 2006; Sakamoto 
2006). GRB 060218 is associated with SN 2006aj and pro- 
vides another example of LL-GRBs. This event is 100 times 
less energetic and the duration is very long ~ 2000 s. A ther- 
mal component in the X-ray and UV-optical spectra was dis- 
covered and the size of the emitting black-body region is es- 
timated to be tbb ~(5x 10 11 - 10 12 ) cm. It would come from 
the shock break out fr om a stellar envelope or a dense wind 
( Cam oana et al, 2006). or from a hot cocoon surrounding the 
GRB ejecta tLiang et al. 2006btlRamirez-ruiz et al. 20 02). 

These GRB 060218-like events are phenomenologic ally pe- 
culiar events compared with the conventional bursts because 
these have lower isotropic luminosity and energy, simpler 
prompt light curves, and larger spectral time lags. Guetta et 
al. (2004) argued that no bright burst z < 0.17 should be ob- 
served by a HETE-like instruments within the next 20 years, 
assuming that GRB 060218-like bursts follow the logN-logP 
relationship of HL-GRBs. Therefore, this unexpected dis- 
covery of GRB 060218 can lead to the idea that they form 
a different new class of GRBs from the conventional HL- 
GRBs, although much uncertainty remains and we do not 
know what distinguishes such LL-GRBs from the conven- 
tional HL-GRBs. Under this assumption, Soderberg et al. 
(2006b) estimated the rate of LL-GRBs and found that they 
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are about ten times more common than conventional HL- 
GRBs, and Liang et al. (2006a) obtained a high LL-GRB rate 
similarly. However, each origin of LL-GRBs could be differ- 
ent. Mazzali et al. (2006) suggested that GRB 980425 and 
GRB 031203 could be related with a black hole formation, 
while GRB 0602 1 8 could be driven by a neutron star. GRB 
0602 1 8-like events might possibly be associated with the birth 
of magnetars. The origin of such LL-GRBs and whether these 
bursts are typical or not are open problems. 

If such LL-GRBs like GRB 060218 are more common than 
HL-GRBs, they would provide enough energy for high en- 
ergy cosmic -rays, neutrinos, and gamma-rays. In this Letter, 
we study the possibility of the high energy cosmic -ray pro- 
duction and the successive neutrino production in LL-GRBs 
under the usual internal shock model. 4 Such neutrino bursts 
from HL-GRBs were predicted in the context of the standard 
scenario of GRBs assuming that ultra-hig h-energy cosmic 
rays (UHECRs) come from GRBs ( Waxman & Bahcal l 19971 
Wax man & Bahcall 1 998). Murase & Nagataki (2006a, 
2006b) also investigated such emission from HL-GRBs and 
from flares, using the M onte Carlo simulation kit GEANT4 
( Agos thTelli et al. 20031) with experimental data. With the 
same method, we study the high energy neutrino emission 
from LL-GRBs. We will also discuss various implications 
for such LL -GRBs. Large neutrino dete ctors such as Ic eCube 
(Ahre nset al. 2004. ANTARE S (lAslanides et al. 1999ft and 
NESTOR (jCirieder et al. 200 ll) are being constructed. In the 
near future, these detectors may detect high energy neutrino 
signals and give us more clues to understanding LL-GRBs and 
testing our model. 

2. THE MODEL 

We suppose GRB 060218-like events as LL-GRBs in this 
Letter. GRB 060218 has low luminosity ~ 10 46 " 47 ergs/s, 
which is much smaller than that of usual HL-GRBs, typically 
£max ~ 10 51-52 ergs/s. Hereafter, we take L max = 10 47 ergs/s 
as a peak luminosity of LL-GRBs and fix E 7 j so = L max 5tN ~ 
jq49-50 ergs as re j ease( j radiation energy. Here 5t is 
the variability time and TV is the number of collisions. To 
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explain the prompt emission, we assume the usual inter- 
nal shock model in which the gamma-rays arise from the 
internal dissipation of ultra-relativistic jets, although there 
is another explanation JPai et al. 2 006). The typical colli- 
sion radius will be expressed by commonly used relation, 
r « 10 15 (r/10) 2 (<5f/150s) cm. Of course, this radius has 
to be smaller than the deceleration radius, r < r^M ~ 4.4 x 

1 /3 

10 16 (£'km.5o/no(r/10) ) cm. The observed light curve of 
GRB 060218 is simple and smooth, suggesting St ~ 10 2-3 s 
dCusumano et al. 2*00 6 ). But it is uncertain whether these pa- 
rameters are typical or not (Fan et al. 2006). Hence, we take 
r ~ 10 14 " 16 cm with T ~ (10-1 00). These radii will be im - 
portant for neutrino production ( Murase & Nagataki 2 006a). 
We also assume that the Lorentz facto r of th e int ernal s hocks 
will be mildly relativistic, T s h ~ (y/Tf/T s + y/T s /T{)/2 ~ 
a few. The typical values in the usual synchrotron model 
are obtained as follows. The minimum Lorentz fac- 
tor of electrons is estimated by 7 c m » e e (m p /m e )(T s h - 
1). Since the intensity of magnetic field is given 

by B = 7.3 x 10 2 Ge^ 1 (r sh (r sh -l)/2) 1/2 4 /2 48 (r/10)- I r7 5 1 , 

the observed break energy is, E h = hj 2 m reB/m e c ~ 
lkeV e 2 e '/ 2 ( r s ,-n 5 / 2 rr. /^i/2,i/2 



= H. 

sh-l) 3/z (r sh /2) l/2 Li I /2 48 r7i, where L M is the 
outflow luminosity. This value is not so different from the 
observed peak energy of GRB 060218, E b ~ keV. 

Although we have too less information about spectral fea- 
tures of LL-GRBs at present, we assume a similar spectral 
shape to that of HL-GRBs for our calculations and approx- 
imate it by the broken power-law instead of exploiting a 
Band spectrum. The photon spectrum in the comoving frame 
is expressed by, dn/de = nb(e/s h ) for e min < e < e° and 

dn jde = n b (e /e b )~^ for e b < e < e max , where we set e min = 0.1 
eV because the sy nchrotron s elf-absorption will be crucial 
below this energy dLi & Song 2004 and e max = 1 MeV be- 
cause the pair a bsorption will be crucial above this energy 
( Asano & Takahara 2003). Corresponding to the observed 
break energy of GRB 060218, E h = 4.9 keV with the assump- 
tion of the relatively low Lorentz factor, we take e b = 0.5 keV 
in the comoving frame as a typical value throughout the Let- 
ter. We also take a = 1 and set (3 = 2.2 as photon indices. Note 
that we may have to wait for other GRB 0602 18 -like events 
to know the reliable typical values. 

We believe not only electrons but also protons will be ac- 
celerated. Although the detail of acceleration mechanisms is 
poorly known, we assume that the first-order Fermi accelera- 
tion mechanism works in GRBs and the distribution of non- 
thermal protons is given by dn p /de p oc e~ 2 . By the condi- 
tion f acc < t p , we can estimate the maximal energy of accel- 
erated protons, where t p is the total cooling time scale given 
by t~ l = f^ + f^yn+fjc+^ad an( ^ tne acceleration time scale is 
given by f acc = r)E p /eBc. Especially, the two time scales f syn 
(synchrotron cooling time) and f a( j ks fd yn (dynamical time) are 
important in our cases. We can estimate the maximum pro- 
ton energy by E p , m . dx « mm[eBr/ r), yj 6we/(T 1 Br](TmpC 2 /m e )] 
from the conditions, f acc < fd yn and f acc < f syn . These two con- 
ditions equivalently lead to, 
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(1) 



These inequalities suggest that the only relatively more lumi- 
nous/magnetized LL-GRBs with higher Lorentz factor (i.e., 
larger Lm and/or e#, and higher T) will have possibilities to 
explain the observed flux of UHECRs. 

We consider neutrinos from the decay of pions generated 
by photomeson productions. The photomeson time scale 
is t pi . Let us evaluate f pl = tdm/tn^ analytically using 
the A-resonance approxi mation (|Murase & Nagataki 2006b; 
IWaxman & Bahcall 1997b as. 



0.06 



v.47 



f {E P /E h p f~ l (E p < E b ) 

r 15 (r/10) 2 £ 5 \ eV \ (E p /E b p r l (E b p <E p ) 



(2) 



where we have used notations such as E p = l0 20 eVE 



p,20- 



where E p ~ 0.5e Ain p c 2 T 2 /E is the proton break energy. 
Here, ea is around 0.3 GeV. From Eq. 0, we can conclude 
that a moderate fraction of high energy accelerated protons 
will be converted into neutrinos. 

Next, we consider the contribution to the neutrino flux from 
a thermal photon component. The discovery of the thermal 
component in GRB 060218 will provide additional photon 
flows. This photon flow has a possibility to produce more 
neutrinos by interaction with protons accelerated in inter- 
nal shocks. We take kT = 0.15 keV and r^s = 10 12 cm as 
the typical photon energ y of the thermal com ponent and the 
apparent emitting radius (Campana et al. 2006), respectively. 
Just for simplicity, we assume the photon density drops as 
oc r~ 2 and approximate it by the isotropic distribution with 
dn/de(e) af«i a b/afeiab(£iabX where dnx&jde^ is the photon 
distribution in the laboratory frame. 

3. RESULTS AND DISCUSSIONS 

We calculate neutrino spectra for some parameter sets and 
will show the case where the width of shells A w (r/2T 2 ) = 
4.5 x 10 12 cm, according to St ~ 150 s. In our calcula- 
tions, we include various cooling processes of pions (syn- 
chrotron cooling, IC cooling, adiabatic cooling) similarly to 
Murase & Nagataki (2006a,2006b) . These cooling processes 
are important for neutrino spectra ( Rachen & Meszaros 1998; 
IWaxman & Bahcall 19971) . A diffuse neutrino background 
under the standard ACDM cosmology (£l m = 0.3, JIa = 
0.7;H() = 71kms _1 Mpc" 1 ) is calculated by using Eq. (15) of 
Murase & Nagataki (2006a), where we set z max = 11. Assum- 
ing that the long GRB rate traces the starformation rate (SFR), 
we exploit the SF2 model of Porciani & Madau (2001) com- 
bined with the normalization of geometrically corrected over- 
all HL-GRB rates ^ H l(0) obtained by Guetta et al. (2005) 
for HL-GRBs. The local LL-GRB rate is very uncertain for 
now. Soderberg et al. (2006b) obtained the geometrically cor- 
rected overall GRB rate, R LL (0) = 230Gpc" 3 yr _1 . Liang et al. 
(2006a) also had a high value, pll(0) = 550Gpc" 3 yr 1 . (Note 
that the true rate /?ll(0) is almost the same as the apparent one 
Pll(0) for LL-GRBs because we are assuming GRB 060218- 
like spherical bursts.) However, too large rates will be impos- 
sible due to constraints by observations of SNe Ibc. Soder- 
berg et al. (2006a) argued that at most ~ 10% of SNe Ibc 
are associated with off-beam LL-GRBs based on their late- 
time radio observations of 68 local SNe Ibc. Hence, the most 
optimistic value allowed from the local SNe Ibc rate will be 
around ~ 4800Gpc~ 3 yr _1 (and the larger value is ruled ou t 
with a confidence level of ~ 90%) (Soderberg et al. 2006a). 
The high rate might be realized if LL-GRBs are related with 
the birth of magnetars and the fraction of SNe Ibc that produce 
magentars is comparable with that of SNe II, i.e., ~ 10%. 
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FIG. 1 . — The neutrino background from GRBs for £ acc = 10 and i 
LL-GRB: r = 9 X 10 14 cm and T = 10 with the local rate ~ 500Gpc" 3 yr' 
obtained by Liang et al. (2006a). LL-GRB (optimistic): r = 9xl0 14 cm 
and r = 10 with the most optimistic local rate ~ 4800Gpc~ 3 yr~' . LL- 
GRB (modest): r = 9 X 10 14 cm and Y =10 with the modest loca l rate 
~ 20Gpc" 3 yr'. HL-GRB: taken from IMurase & Nagataki 20063) with 
E^. isn / N = 2 X 10 51 ergs, r= 1Q13-1 4.5 cm p = 300 wg. Waxman-Bahcall 
bounds iWaxman & Bahcall 1998). £b and § acc are the ratio of energy den- 
sity, £b = Ub/Uj and £ acc = U p /Uj, respectively. For the fast cooling case 
and the acceleration efficiency ~ 1, we have £g ~ and £ acc ~ l/e f . 



Although we calculate numerically, we can estimate the 
diffuse neutrino flux from LL- GRBs a pproxima tely by the 
following analytical expr ession (IMurase & Nagataki 20 06b: 
IWaxman & Bahcal l 19981) . 



4ttH 4 
7 x 10" 10 ' 

0.05 



1 mm[\J pl ]E 2 p ^R^{0)f z 



GeVcm-V'str- 1 feW 



50 



^ll(O) 



500Gpc~ 3 yr-! 
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where £ll is the geometrically corrected radiated energy of 
LL-GRBs, f z is the correction factor for the possible contribu- 
tion from high redshift sources, and we have used e 



p. max 



10 v 

GeV. Our numerical results are shown in Fig. 1. From these 
results, we can estimate the number of muon events due 
to muon-neutrinos above TeV energy by using Eq. (18) of 
Ioka et al. (2005) as the detection probability and a geo- 
metrical detector area of Ad et = 1km . From Fig. 1, we 
can obtain = 6.6 events/yr for Pll(0) = 500Gpc~ 3 yr _I and 
N/j, = 64 events/yr for the most optimistic local rate. We also 
show the modest case where the local rate of LL-GRBs is 
comparable to the geometrically corrected local rate of HL- 
GRBs, Rhl(0) ~ 20Gpc" 3 yr 1 . In this case, we can find 
N/j, = 0.3 events/yr. The neutrino backgrounds from LL-GRBs 
have possibilities t o be com parable with that f rom HL-GRBs, 
N fl = 17 events/vr ( Mur ase & Naga taki 20063). 

Unfortunately, neutrino signals from LL-GRBs are dark in 
the sense that most signals will not correlate with the prompt 
emission. Only for very nearby bursts, we might be able to 
expect their correlations and it will need many-years opera- 
tions. The BAT detector on Swift has the sensitivity to de- 
tect the bursts > 10~ 8 ergscm" 2 s _1 . Hence, we can expect 
correlated events only when < 300 Mpc for bursts with 
imax ~ 10 47 ergs/s. The expected correlated muon events are 



~ 1 events per 1 1 years for Pll(0) = 500Gpc~ 3 yr _1 . 

However, SNe Ibc associated with LL-GRBs could be de- 
tected by optical-infrared follow-ups triggered by a neutrino 
event. The angular resolution of IceCube for neutrinos is 
about 1 degree or so, which can be searched with wide- 
field cameras such as Suprime-Cam on the Subaru telescope 
(whose field-of-view is 0.5 degrees) up to z ~ 1.2. In the 
field-of-view we would find ~ 10 SNe and ~ 1 SNe Ibc that 
exploded within ~ 1 month. With the SN light curves ~ 10 
days after the burst, we can pin down the burst time within ~ 1 
day or so, during which the atmospheric neutrino background 
within 1 degree would be small, i.e., < 0.1 events/day for 
abo ve TeV energy neu trinos and less for higher energy thresh- 
old jAndo & Beacom 20"05l) . In addition, SNe Ibc could be 
specified by using telescopes such as HST Therefore, we can 
in principle detect LL-GRB neutrino events associated with 
SNe Ibc, even though X/7-rays are not observed by Swift. 
The expected number of muon events is =2.4 events/yr 
for LL-GRBs within z ~ 1.2, with r = 9 x 10 14 cm, T = 10, 
and p LL (0) = 500Gpc- 3 yr _1 - O f course, such a follow-up with 
SNe detections will be difficult and it is severer to distinguish 
SNe Ibc from SNe la at higher redshift. Nevertheless, it is 
worthwhile to develop this kind of possibility of high energy 
neutrino astronomy not only for finding far SNe Ibc associ- 
ated with LL-GRBs but also for revealing their origins. 

We can expect high energy neutrinos from one LL-GRB 
only if the burst is nearby or energetic, similarly to the case 
of HL-GRBs. In Fig. 2, we show an example of the observed 
neutrino spectra from the source at 10 Mpc. The expected 
muon events from neutrinos above TeV energy are = 1.1 
events in the case of T = 10 in Fig. 2. If we can detect such an 
event, we will be able to obtain some information on £ acc , £b, 
the photon density, the duration of bursts, and so on. In Fig. 
2, we also show the contribution from the thermal target pho- 
ton. The GRB 060218-like bursts could provide us = 0.2 
events originating from the interaction between nonthermal 
protons and the thermal photon flow. This result depends on 
the temperature of the black body region. 
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FIG. 2. — The observed muon-neutrino (v^ + D^) spectra for one very 
nearby GRB event at 10 Mpc. Solid line: r = 9 X 10 14 cm and V = 10. Dashed 
line: the contribution from the blackbody target photon with r = 9 X 10 14 cm 
and T = 10. Dotted line: r = 2.25 X 10 14 cm and Y = 5. Dot-dashed line: 
r = 3.6 X 10 15 cm and F = 20. In all cases £g = 1 and £ acc = 10 (see the 
caption of Fig. 1). Note, the case of T = 20 would not be plausible because 
the magnetic field strength seems too small to explain the prompt emission 
by the standard model. 
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HL-GRBs may be the main sources of UHECRs 
iWaxman 19951) . The optical thickness for the photomeson 
production can be smaller than unity especially at larger radii 

1 /2 

r > 10 14 CE 7 ,i so /M0 51 ergs) cm in the internal shocks of 
HL-GRBs and UHECRs can be produced in such regions. 
In the case of LL-GRBs, it seems more difficult to accel- 
erate protons up to ultra-high energy due to the lowness of 
their luminosities and it would need some fine tuning (see Eq. 

although we have to know about their properties such as 
their luminosity function. Even if the acceleration to ~ 10 20 
eV is difficult, the energy budget of LL-GRBs could be large 
enough to explain UHECRs (~ 10 44 ergsMpc~ 3 yr _I ) because 
the cosmic -ray production rate per Mpc 3 is estimated by, 



,dN B 



E p 7^ ~ 2 - 5 x 1 ° 43 er § s m p c " j y r " 

at p 



1 f £ac 



10 



x NL m!ixA7 r 15 [ — 



Pll(0) 



SOOGpc^yr" 1 



(4) 



18.5 



Therefore, when the maximum proton energy exceeds 10 
eV, the neutrino flux should be constrained by the observed 
flux of UHECRs, even if LL-GRBs cannot explain all UHE- 
CRs. This implies that we have possibilities to constrain phys- 
ical parameters of LL-GRBs by the observation of UHECRs. 

High energy neutrino emission cannot avoid high energy 
gamma-ray emission through the neutral pion decay. Such 
high gamma-rays would cascade in the source and/or in mi- 
crowave and infrared background (iDermer & Atovan 20041 
iRazzaaue et al. 20041) . The detailed calculation is needed to 
calculate the expected spectra. However, the detection of such 
high energy emission by GLAST and/or the BAT detector on 



Swift would be difficult except for nearby and/or energetic 
bursts similarly to the cases of neutrinos. 

If there is a shock break out which might be the origin of 
the mysterious thermal component, the shock may become 
collisionless, and protons may be accelerated there as well 
as electrons, so that neutrin os could be produced through 
the pp interaction ( Wax man & Loeb 20011) . In addition, pro- 
tons accelerated by internal shocks inside the stellar envelope 
could prod uce detectable ~ TeV neutrinos mainly by pp in- 
teracti ons (Meszaros & Wax man" 200 It IRazzaaue etal. 20031: 
lAndo & Beacom 2005). These other neutrino signals might 
become clues on the connection between GRBs and SNe. 

In this Letter, we have discussed a possibility that LL-GRBs 
could produce UHECRs and detectable high energy neutri- 
nos. Because the possible higher rate of LL-GRBs can cover 
the relatively lower energy of them, we can expect that the 
diffuse neutrino flux could be comparable with that of HL- 
GRBs. Of course, the results depend on several unknown pa- 
rameters such as the bulk Lorentz factor, and future observa- 
tions are needed for more realistic predictions. If parameters 
we have adopted are typical, expected neutrino signals may 
not only give us independent information but could also be 
useful as one of indicators of far SNe Ibc and may contribute 
to revealing their mysteri ous origin. One of possibilities is the 
birth of magnetars (Soderberg et al. 2006b). Possibly, metal- 
licity might play a crucial role ( Stanek et al. 2006). 

We thank the referee, M. Doi, N. Tominaga, R. Romani, 
S. Jha, K. Toma, N. Kawanaka, J.F. Beacom, and E. Liang 
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for Scientific Research of the Japanese Ministry of Education, 
Culture, Sports, Science, and Technology 18740147 (K.I.). 
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